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We establish that the entropy production rate of a classicaUy chaotic Hamiltonian system coupled 
to the environment settles, after a transient, to a meta-stable value given by the sum of positive 
generalized Lyapunov exponents. A meta-stable steady state is generated in this process. This 
behavior also occurs in quantum systems close to the classical limit where it leads to the restoration 
of quantum-classical correspondence in chaotic systems coupled to the environment. 
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0^ , The Lyapunov exponents of a chaotic classical Hamiltonian system characterize the rate of exponential divergence 

■ of neighbouring phase-space trajectories thus defining the microscopic time-scales. The question [Q of their role in the 
I : macroscopic behavior of the system is hence at the foundation of statistical mechanics. It is also of deep interest to 

understand how the Lyapunov exponents govern the behavior of the quantum counterparts of these chaotic systems 
1^. For instance, it is predicted ||] that quantum effects comparable to the classical behavior should be measurable 
on a time-scale ts ~ J^^iji^- Here h is Planck's constant scaled by a characteristic action for the system and A 
is the largest Lyapunov exponent of the system. This 'break-time' can be w 20 years for astronomical systems 
' and is clearly incompatible with observations. It is argued |^ that this problem may be resolved by considering the 
' , decoherence resulting from including environmental effects on the system evolution. Since this can be understood as 
^ ' a dynamical coarse-graining, and statistical properties of the system are typically obtained through a static coarse- 
^ , graining, a diagnostic of common primary interest for both issues is the behavior of the coarse-grained entropy of the 

■ system. The Gibbs entropy Sg = Tr[plnp] and the linear or Renyi entropy ^2 = ln(Tr[p^]) are both useful j|] in this 
. regard. Here p is the classical phase-space probability distribution {x) or the quantum quasi- probability Wigner 

■ distribution p^ {x) as appropriate and Tr denotes the integration over all phase-space variables x — {q,p). These 
^ I , entropies measure the information in the probability distribution and remain constant for Hamiltonian evolution 

■ without coarse-graining. A change in the coarse-grained entropy corresponds to information being lost at the finest 
, scales. Since in a chaotic system p acquires structure at small scales exponentially rapidly ||^, S2 typically grows 

■ exponentially rapidly in magnitude initially for such classical systems, as also for their quantum analogs close to the 
classical limit 0. Numerical studies with static coarse-graining |^ and the analysis of the behavior of a Gaussian 

'■^ ■ wave-packet in the upside-down harmonic oscillator suggest that for weakly coupled systems Sg settles to a meta- 
Q : stable 'constant' after this transient. This meta-stable value is conjectured [Q to be the sum over the positive 

' Lyapunov exponents of the classical system, independent of the details of the original distribution. This conjecture 
f~| , has not yet been verified. Preliminary analyses ||] indicate that systems may become meta-stable although the 
O ' quantitative conjecture does not hold in general. In this paper, we consider a variation of this conjecture and are able 
^ , to analytically establish it for the general case of an arbitrary distribution evolving in a chaotic system. We do so by 
considering the dynamics of the quantity x{t) which measures the degree of structure in a phase-space distribution, 
and is directly proportional to the entropy production rate. We demonstrate that under reasonable approximations 
• the entropy production rate S2 indeed settles after an exponential transient to S2 ~ ~^i-^2i- ^2,i, as defined 
below, are generalized Lyapunov exponents that depend upon the initial distribution in general. We further show 
that this is a robust result: Although for sufficiently weakly coupled systems 5*2 increases in magnitude to a constant, 
for a stronger coupling the magnitude of S2 starts high and decreases to the same constant. The meta-stable rate is 
hence independent of the environmental strength. We also generalize the idea that for such systems p can evolve to a 
meta-stable state and show that for such states the quantum corrections for classically chaotic systems remain 

small for all times, indeed resolving the issue of the anomalously short 'break-time'. We thus obtain a condition for 
the restoration of correspondence through decoherence. Finally, we present numerical results verifying these analytic 
results for the particular case of the quantum Cat Map and its classical limit, a uniformly hyperbolic K-system 
satisfying the assumptions in our arguments . 

The evolution of p^ under the potential V{q) and coupled to an external environment is 



9^ ^ SH 4- V o^"+V^ , ^^2w 



The first term on the right is the Poisson bracket, generating the classical evolution for p^; the terms in h add 
the quantal evolution. The environmental effects are contained in the term; for simplicity we have assumed 
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coupling to all the phase-space variables, although the results can be easily generalized. The parameter D depends 
strongly upon the form of the coupling and the spectrum of the environment [T^ . The entropy production rate is 

then S2 = 2D ^^^^^^^^ ' — ~2D j,l.'^^^wyi^ = — 2I?x where the second equality results from an integration by parts. 
The global quantity x^(^) corresponds to the mean-square radius of the Fourier expansion of p; it hence measures 
the structure in the distribution and also governs the entropy production rate. To understand the behavior of 
S2 first consider the classical Hamiltonian evolution, that is, only the Poisson bracket terms in Eq. ([|). In this 
case, the point dynamics satisfy x — f{x) with the accompanying induced flow xt — $*xo and the phase-space 
gradient along a trajectory is governed by p{xt) = — A4(t, a;o)Vp(a;o), where the fundamental matrix Ai is given 
by |l[ the time-ordered series A^(t, Iq) = Texp J^* ^dr. This defines a real symmetric matrix Ai"^ ■ M (the 

transpose denoted by M^) which can be diagonalized as M'^{t,xo) ■ M{t,xo) = '^^Ui{t,Xo)ai{t,xo)uJ {t,xo), 
whence the Ui constitute a local orthonormal tangent space for the flow. The local Lyapunov exponents are then 
given by Ai(a;o) = limt^oo jj- In cri(t, xq). We may now define the global averages 



A2,j(p,t) = ^In 



Tr[|V,p(xo)|2a,(t,xo)] 



Tr[V,p(xo)|2 
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Tr[|V,p(a;t)| 
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Tr[V,p(xo)|2 



(2) 



and their limits A2,i(p) — limt-^oo ^2a{p, t). The last equality in Eq. (j^) comes from the orthonormality of the m and 
the dependence of Vp(a;t) on A4. These A are p dependent versions of the usual generalized Lyapunov exponents 
of second order |^ which may be recovered by replacing \\7ip{xo)\'^ by the natural invariant measure po{x) of the 
dynamics. For linear systems ai(t,xo) = cri(t) and the definitions are independent of p. For ergodic dynamics 
^2,iip,t) = A2,i(/5) for all time. However, numerical evidence suggests that this equality is valid within small errors 
for 'typical' initial densities p for other Hamiltonian flows as well. We can always write A2^i{p,t) — K2.i{p) + Ci/t 
where ^ fluctuates in general and ^/t vanishes with increasing t. The analysis so far has been exact; hereafter we make 
the approximation of neglecting the ^ term. We now decompose the entropy production rate into the contribution 
from the different stability subspaces as ^2 = -2I?x^ = ^2_D^- xf where xf — Tr[|Vip(a;t)p]/Tr[p^(a;f)] = {^i)p'l 
the i component of x^. The preceeding enables us to write 

^ « 2x?A2, (3) 

or equivalently 5*2 « — 2L) xf (0)e^^^''* where the p-dependence of A2,i is left implicit hereafter. This shows the 
initial behavior of 5*2; it is dominated by the largest positive generalized exponent A^j^, thus recovering an earlier 
approximate result [Q. After this initial stage, the structure reaches the finest scales and hence the D-dependent 
diffusion term becomes important. To include the impact of this we first consider purely diffusive behavior such that 
dtp — D'S/^p. Since the Laplacian is independent of the coordinate system, this may be decomposed into the same 
subspaces as above as dtp = Dj^i ^iP with the solution being a product of the distributions as p = IliPi{xi). Each 
direction has the standard solution written in terms of the Fourier components as p^. = pk^ exp(— £)fc|i) whence the 
individual xK*) = i^i)pl = Z^fc exp(-2_Dfc^t)/(^j. |pfc(0)|^ exp(-2_Dfc^i)) have the time-dependence 

^ = ^2D[{kf)^. - {kf)l.J « -4D(fc?>2, = -4Dxt. (4) 

The approximation in Eq. (Q) is a mean-field one, valid for the usual Gaussian solution to the diffusion equation for 
example. The various xf therefore behave as follows: An initial exponential transient entirely kills these quantities 
in the stable directions (corresponding to the negative Lyapunov exponents) and since this is enhanced by diffusive 
effects, these directions have a negligible role thereafter. In the unstable directions, however, the initial exponential 
growth is balanced by the diffusion. If we now explicitly set dxi/dt equal to the sum of the chaotic [Eq. ^ and 

diffusive [Eq. ^ terms, we get = 2A^jxf ~ 4Dxt for the unstable directions. This has the stationary solution 

X^* = A^-/2D yielding the stationary entropy production rate ^2* = — 2i?x^* = ~2D^jX?* — ~J2i^2i '^ith the 
sum over the positive exponents only, as just argued. Thus, within the approximations as above, we have shown that 
the entropy production rate for a chaotic system weakly coupled to the environment settles after an initial exponential 
transient to a meta-stable value as above, independent of the precise magnitude of the environmental effects [ 1^ . This 
solution is stable under small pertubations; thus, the exact dynamics arguably leave it unaltered. We now consider 
possible constraints on this result. First, note that the transition from the exponential to the linear regimes happens 
at t* « (l/2A^j) ln[A^j/2Z?x|(0)] [the shortest scale is approximated by the largest exponent A^j^ = A and the full 
X^ and we shall use those hereafter]. Since = for /3 = 2A~^I?x^(0) = 1, for 3> 1, the initial behavior of Xi is not 
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the exponential transient but is given by Eq. instead. For this case, the initial diffusive effects are balanced by the 
chaos such that an initially large 5*2 decreases to a constant. Second, the difference between the initial entropy 52(0) 
and the final entropy 6*2(00) is typically finite. Thus, the system may not have 'enough initial entropy' to evolve to 
the constant entropy rate as above. A rough estimate (and numerical results as described below) indicate that we 
need 82(0) — 62(00) > 3 for the linear saturation behavior to emerge iQ. Finally, in contrast to this behavior for 
chaotic systems, a similar analysis for non-chaotic systems yields an initial non-exponential decay of 62 after which 
the competition from the diffusive term implies a maximum (oc D~^) for 62, and a rapid decrease thereafter. As 
a result of these dynamics, appropriate initial distributions p for chaotic systems settle to a meta-stable state with 
as above. This is a remarkable result, generalizing earlier arguments for a 'steady-state Gaussian width'. 

Further, this enables us to bound the quantum corrections in Eq. (|^). Consider the first quantum term in Eq. (|^), 
which scales as Ti^x^- In the absence of the Z?-dependent terms, x grows exponentially rapidly as exp(At), thus leading 
to a quantum-classical 'break-time' w i ^^(l) after which the quantum 'correction' is comparable to the classical 
evolution and a classical description of the dynamics is invalid. This logarithmic dependence of on ti is extremely 
weak and is a point of debate in the analysis of quantum chaotic systems jj] . However, environmental effects saturate x 
and hence the first quantum term at C ~ h^{X/2D)^/^ such that this indeed is a 'correction' of Oih^) to the classical 
evolution and there is no 'break-time'. The condition [p^ for restoration of quantum-classical correspondence for 
chaotic systems may be summarized as C <C 1 (note that the other quantum terms are higher order in C) . Physically, 
quantum effects are important at the smallest scales k T% oi phase-space. Classically chaotic evolution increases the 
support of a distribution at the finest scales exponentially rapidly, thus enhancing quantum effects. However, noise 
or coarse-graining washes out the details of the fine-scale structure, thus restoring quantum-classical correspondence 
For systems where the above inequality is violated, however, a classical or semiclassical analysis breaks down 
and we must use the exact quantum evolution. 

We have numerically verified the theory using the equivalent of Eq. ([l]) for the classical and quantum Cat Map . 
As this system is linear and ergodic, the approximations above hold exactly. Some of the data obtained are shown 
in Figs. (1-3). First, with a scaled Planck's constant 2-111% = 10~^, we expect essentially classical behavior, even 
with no added noise. This is indeed verified as in Fig. (1) where the classical data for the same initial conditions 
and D exactly overlie the quantum results. For an initially sharply localized Gaussian state, corresponding to high 
initial entropy, when D is small we see initial exponential entropy production |^]; this saturates to linear behavior 
with the predicted slope. The initial exponential behavior is swamped as D is increased, as argued previously. For 
a second initial state which is spread out and has lower initial entropy, we see similar stable behavior - although 
62 now saturates earlier, as discussed above - confirming that this meta-stability is essentially independent of initial 
conditions and of D for near-classical systems. Fig. (2) shows the quantum-to-classical transition by comparing 
various runs for a single sharply localized initial condition at 27r?i = 10^^ with their classical equivalents. We see 
that for very small D (hence large C) there is significant quantal deviation from the classical behavior. However, as 
D is increased (and C correspondingly decreases below 1), the quantal 62 approaches the classical linear saturation 
behavior as above. A more spread-out initial distribution has the same qualitative behavior (not shown) except 
that it saturates early as expected. We note that the distributions for both Figs. (1,2) are not Gaussian during the 
meta-stable stage of the dynamics. Fig. (3) shows results for 27r?i = 0.1. Here the uncertainty principle constrains 
the initial p to occupy a substantial portion of the available phase-space. We see again that for C 1, the quantum 
behavior is substantially different from the classical and there is almost no sensitivity to the environment (low entropy 
production). As D increases, quantum-classical correspondence is indeed restored, although D is now so large the 
dynamics are essentially that of the noise alone and 62 saturates at less than A2. The abruptness of the tail in all 
our computations is most probably a numerical artifact due to the extremely small numbers being computed at that 
time; a detailed understanding of this tail is still absent. These results show that the entropy production rate for an 
arbitrary distribution in classically chaotic system saturates, after an exponentially rapid transient, to the sum of the 
positive generalized Lyapunov exponents of the system; further, the distribution settles to a meta-stable state. This 
is a clear signature of the underlying chaos, stable against perturbation. This behavior is echoed by quantum systems 
close to the classical limit, and the saturation is the precursor for quantum-classical correspondence in chaotic systems 
coupled to the environment. 
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FIG. 1. Time-dependence of the entropy 5*2 (f) for the chaotic Cat Map for varying levels of coupling to the environment, as 
measured by D. Classical and quantal {211% — 10"'') numerics were run for identical initial conditions and noise strengths D. 
The classical results are shown as points overlaid on the lines for the corresponding quantal data. There are two different initial 
conditions, one with high initial entropy ^2(0) and the other with lower ^2(0). A reference line with slope equal to A2 — 0.9624 
is also shown. 

FIG. 2. Quantal and classical results are shown with 27vh = 10"^, a single initial condition and varying D. As D is increased, 
the quantal behavior approaches that of the classical system. As in Fig. (1), a reference line with slope equal to A2 = 0.9624 
is also shown. 

FIG. 3. As in Fig. (2), but with 2Trh = lO"'^. The quantal system is almost unaffected by small D. The quantal behavior 
tends to the classical with increasing D; however, it is then dominated by noise and the slope is always less than that of the 
reference line shown with slope of A2 = 0.9624. 
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